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1 3.  ABSTRACT  {Maximum  200  vjords) 


A  jet  diffusion  flame  attached  to  a  burner  rim  may  lift  off  and  become  stabilized  further  downstream  when  the  jet  velocity  is 
sufficiently  increased.  In  turbulent  jet  diffusion  flames,  such  liftoff  has  been  described  alternately  as  the  result  of  a  stabilized  premixed 
flame  base,  and  as  the  result  of  extinguished  diffusion  flamelets  at  the  flame  base.  Laminar  flames  exhibit  liftoff  behavior  as  well,  and 
possess  a  relatively  simpler  flame  structure  which  may  be  studied  to  provide  insight  into  the  basic  mechanism  responsible  for  flame 
lift.  In  the  present  study,  an  asymptotic  solution  and  a  numerical  solution  of  a  reduced  set  of  equations  are  used  to  study  the  lifted 
structure  of  a  laminar  diffusion  flame  associated  with  a  fuel  jet.  The  numerical  model  solves  a  temperature  equation  with  a  finite 
chemical  reaction  rate.  The  asymptotic  analysis  is  based  on  a  flame-sheet  model  in  the  limit  of  large  activation  energy  and  large 
Damkdhler  number.  While  it  is  likely  that  premixing  at  the  flame  base  affects  the  structure  of  a  lifted  flame,  the  current  analysis 
suggests  that  local  extinction  of  the  flame  at  its  base  is  a  contributing  mechanism  to  flame  lift. 
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AN  ANALYSIS  OF  LIFT-OFF  IN  LAMINAR  DIFFUSION  FLAMES 


INTRODUCTION 

When  the  flow  rate  of  a  fuel  jet  is  increased,  a  jet  diffusion  flame  attached  to  a  burner  rim 
may  lift  oflf  and  become  stabilized  further  downstream.  Theories  describing  the  stabilization  of 
such  lifted  flames  in  turbulent  flow  regimes  have  postulated  that  the  flame  at  its  base  is  fully 
premixed  and  thus  stabilized  where  the  average  axial  velocity  equals  the  mean  burning  velocity  of 
the  premixture  [1],  [2].  Other  theories  suggest  that  turbulent  diffusion  flames  may  be  described 
as  a  collection  of  laminar  flamelets,  with  liftoff  occurring  when  a  sufficient  number  of  flamelets 
near  the  jet  exit  become  extinguished  due  to  excessive  strain  [3],  [4].  An  extensive  review  of  such 
turbulent  liftoff  theories  has  been  given  by  Pitts  [5].  Laminar  flames  exhibit  liftoff  behavior  as 
well,  and  have  a  relatively  simpler  flame  structure  which  may  be  studied  to  provide  insight  into 
the  basic  mechanism  responsible  for  flame  lift.  In  the  particular  context  of  a  fuel  jet,  lifted  laminar 
flames  have  been  observed  both  experimentally  and  numerically  [6],  [7].  Theoretical  predictions 
of  laminar  flame  liftoff  heights  have  been  obtained  using  cold  jet  theories  in  which  the  flame  base 
is  predicted  to  occur  along  stoichiometric  contours  where  the  flow  velocity  balances  the  velocity 
of  a  premixed  flame  base  [8],  [6],  [9].  Wichman  and  Ramadan  have  used  scaling  arguments,  and 
detailed  comparisons  to  numerical  solutions  of  the  same  equations,  to  elucidate  liftoff  trends  [10]. 
Experimental  and  numerical  studies  which  may  be  relevant  to  the  flame  structure  at  the  base  of  a 
lifted  flame  include  those  of  triple  flames,  often  studied  theoretically  in  the  context  of  a  stratified 
premixture  [11],  as  well  as  those  describing  diffusion  flame  evolution  in  mixing  layers  [12]  and 
phenomenological  models  of  edge  flames  [13]. 

In  the  present  study,  the  structure  of  a  laminar  diffusion  flame  associated  with  a  fuel  jet  is  ana¬ 
lyzed  to  investigate  flame  lift.  Here  a  Cartesian  geometry  is  considered  in  which  the  two-dimensional 
planar  fuel  jet  is  surrounded  by  a  coflow  of  oxidant  of  the  same  velocity.  Two  approaches  are  under¬ 
taken  to  describe  such  a  flame.  First,  solutions  are  found  by  solving  a  reduced  system  numerically. 
Next,  an  asymptotic  solution  is  constructed.  The  formulation  of  the  asymptotic  model  is  based 
on  a  flame  sheet  approximation,  in  a  distinguished  limit  between  large  activation  energy  and  large 
Damkohler  number.  Results  of  the  two  solution  methods  are  compared  and  discussed.  The  effect 
of  flow  velocity  as  well  as  other  parameters  of  the  combustion  field  on  the  flame  is  investigated. 

FORMULATION 

The  jet  diffusion  flame  is  modeled  by  a  two-dimensional  slot  geometry,  as  shown  in  Figure  1. 
The  half-width  of  the  inner  wall  has  been  normalized  to  unity.  In  order  to  simulate  an  unconfined 
flame,  the  half  width  of  the  (normalized)  outer  boundary  L  is  taken  to  be  large  so  that  the  outer 
wall  has  a  negligible  effect  on  the  combustion  field.  The  chemical  reaction  is  modeled  as  a  global 
one-step  reaction  of  the  form 

vpF->rvxO-^P  (1) 

where  F,  O  and  P  represent  the  fuel,  oxidant  and  product  of  the  reaction,  and  vpiVx  are  the  fuel 
and  oxidant  stoichiometric  coefficients.  An  Arrhenius  temperature  dependence  of  the  reaction  has 
been  assumed. 
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The  resulting  dimensionless,  steady-state  form  of  the  governing  equations  reflecting  conservation 
of  energy  and  mass  balance  of  the  reactants  is 


dT  d^T  d^T  , 

(2) 

_dYp  d^Yp  d^Yp  ^  2,, 

”  dx 

(3) 

dYo  d^Yo  d^Yo  ^  2,,  ,,  «/t 

”  to  -  ax^  -  -  -Dp^YoYce-‘IP 

(4) 

Here  it  has  been  assumed  that  the  velocity  profile  is  one-dimensional, 

V  =  («,0,0) 

(5) 

so  that,  from  the  continuity  equation,  the  mass  flux  m  =  pw=constant.  The  momentum  equation 
is  thus  decoupled  from  the  system  and  can  be  solved  a  posteriori  to  determine  the  small  deviation 
from  the  constant,  time-independent  ambient  pressure.  The  equation  of  state  is  given  by  pT  —  1. 

In  the  above  equations,  lengths  have  been  made  dimensionless  by  the  half  width  of  the  fuel 
slot  a,  and  velocity  by  the  diffusional  velocity  =  \j p^c-pa.  Here  Cp  is  an  averaged  constant 
specific  heat  of  the  mixture,  and  A  is  the  thermal  conductivity.  A  reference  value  pc  =  PcW/RqTc 
is  used  to  make  the  density  dimensionless,  where  Pc  is  the  ambient  pressure,  W  is  an  averaged 
molecular  weight  of  the  mixture,  and  Rq  is  the  gas  constant.  The  temperature  T  has  been  scaled 
=  {Q/cpjiyFj/t^FWp),  the  fuel  mass  fraction  by  its  initial  value  at  the  slot  exit  Ypj,  and 
the  oxidant  mass  fraction  has  been  scaled  by  Ypj{uxWx /FpWp).  Here  Q  is  the  heat  release,  and 
Wp  and  Wx  are  the  fuel  and  oxidant  molecular  weights,  respectively.  Unity  Lewis  numbers  have 
been  assumed,  so  that  A/ pCpDp  =  XJ pCpDx  =  1  where  Dp  and  Dx  are  fuel  and  oxidant  mass 
diflFusivities.  Then  the  Damkohler  number  D  =  FxYpjBpco/Wpu^),  where  5  is  a  pre-exponential 
factor  of  the  reaction,  and  the  scaled  activation  energy  0  =  E/RqTc- 

The  configuration  is  symmetric  with  respect  to  y,  so  that  it  suffices  to  solve  the  system  on  the 
half-space  0  <  y  <  L.  Boundary  conditions  reflecting  symmetry  at  y  =  0  and  no-flux  conditions  at 
y  =  L  are 

dT  dYp  dYo 

dy  dy  dy  ~ 

At  a;  =  0,  boundary  conditions  are 

=  0.  for  0  <  y  <  1  (7) 

V  ^yo  ^  rr.  dT  dYp 

''^Yo  —  mXoo,  rnT  -  —  =  mToo,  rnYp  — =  0  for  1  <  y  <  £,,  (8) 

where  Tj  and  Too  are  the  (dimensionless)  temperatures  of  the  fuel  and  oxidant  at  their  sources, 
respectively.  The  dimensionless  oxidant  source  mass  fraction  is  Xqo  and  the  scaled  fuel  source  mass 
fraction  is  unity.  Boundedness  of  solutions  is  required  as  a;  — >■  oo. 

Before  proceeding  to  finding  solutions  of  the  above  equations,  we  note  the  existence  of  coupling 
functions  which  satisfy  a  reaction-free  form  of  the  governing  equations.  In  particular, 

L{T  +  Yp)  =  L{T  +  Yo)  =  LiYp-Yo)  =  0  (9) 

where 

^  _  d  d^  d^ 

^-"^dx  dx^ 
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Applying  the  boundary  conditions  of  Eqs.  6-8  one  thus  finds  solutions 

'J'  _L  ^  "f"  '^Oo{L  ~  1)  .  r.  — 

^  ~  L  +  [1 +  -  7bo]5(a;,y) 

ra.v„  _  Tj+{Too  +  Xoo)(L-\) 

- - - L~^ - -  +  [Tj  -  Too  -  Xoo]S(x,  y) 


Yf-Yo  =  ^^^^^^^^-^  +  [l  +  Xoo]S{x,y) 


where 
S{x,y)  = 


(11) 

(12) 

(13) 


4m 


, _ _ _ _  ( nF\  ( nwy 

nn{m+  +  A[nF/LY)  V  X  j  (.T 


/ 

m  calculating  S{x,y),  we  truncate  the  infinite  series  to  a  finite  series  over  N  modes  The  above 
solutions  will  prove  useful  in  simplifying  the  analysis  which  follows. 

NUMERICAL  FORMULATION 

Pn  Jo  w'"  numerically  solve  a  time-dependent  form  of  the  system 

yI  and  y  "•  to  a  steady  state.  The  coupling  functions  Eqs.  11-12  may  be  used  to  express 

equation  °  °f  ‘h'  temperLre 

dT  dT  d^T  d^T 

=  Dp^Y^yoe-‘i^, 

subject  to  the  boundary  condition  at  ar  =  0 


(14) 


forO<t/<l,  mT-^  =  mTj;  for  1  <  y  <  L,  mT-~  =  mToQ 


and  the  condition,  at  a;  =  3:^02;, 
At  y  =  0  and  y  =  L, 


dT 

dx  “°- 


dy 


=  0. 


(15) 

(16) 
(17) 


Equation  14  is  solved  using  a  flux-corrected  transport  method.  The  liftoff  height  is  calculated 

7’iiftoff  =  (ig) 

temUJturrs*'7\J'”mlxTJ%^^^^  maximum  of  the  fuel  and  oxidant  source 

buXn  Tl!!’ J"™  r  ?>- ^00);  Tj  IS  the  flame  temperature  corresponding  to  complete  com¬ 
bustion.  The  flame  height  in  the  numerical  calculations  is  defined  to  be  the  x-location  where  Je 
temperature  reaches  a  maximum  value  at  the  centerline  y  =  0.  Converrence  to  a 
determined  primarily  from  the  convergence  of  the  flame  liftoff  height.  Convergence  tests  were  all^ 
performed  with  respect  to  grid  discretization  and  timestep  size. 
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Numerical  Solutions 

Typical  profiles  resulting  from  the  solution  of  Eqs.  14-17  are  shown  in  Figure  2.  Here  Tj  =  Too  =  0.1, 
Xoo  =  1.0,  Z)  =  4  X  10',  and  9  —  5.0.  The  liftoff  height  in  Figure  2  is  calculated  to  be  x^h  =  2.85, 
and  the  flame  height  is  Xfh  =  9.93. 

In  Figure  3,  liftoff  heights  versus  flow  rate  m  are  shown  for  several  values  of  the  Damkohler 
number  D.  Here  Tj  =  Too  =  0.1,  Xoo  =  1-0,  and  $  =  5.0.  In  this  figure,  the  liftoff  height 
exhibits  a  threshold  phenomenon  with  respect  to  the  flow  rate  parameter  m.  The  flame  appears  to 
remain  attached  until  a  sufl[iciently  large  value  of  m  is  reached;  as  m  is  increased  beyond  this  value 
the  liftoff  height  increases  at  a  faster-than-linear  rate.  Calculated  flame  heights  are  also  shown  in 
Figure  3,  and  increase  linearly  with  flow  rate  m. 

In  Figure  4,  liftoff  heights  are  shown  for  several  values  of  m  as  a  function  of  the  Damkohler 
number.  Here,  as  in  Figure  3,  Tj  =  Too  —  0.1,  Xoo  =  1-0,  and  9  =  5.0.  Figure  4  again  indicates 
the  trend  that  as  the  Damkohler  number  decreases,  liftoff  heights  increase  for  a  fixed  value  of  to. 
Furthermore,  at  a  given  value  of  D,  liftoff  heights  are  higher  for  larger  values  of  to. 

In  Figure  5,  temperature  profiles  are  shown  for  three  values  of  the  flow  rate  to.  Here  Tj  =  Too  =  0.1, 
Xoo  =  1-0,  D  =  4x10^,  9  =  5.0.  In  Figure  6,  corresponding  maps  of  the  reaction  rate  p^DYFYoe~^l^ 
are  shown.  As  to  is  increased,  the  flame  evolves  from  a  Burke-Schumann  type  attached  flame,  to 
a  lifted  flame  with  triple  flame  structures  at  its  base,  to  a  higher,  lifted,  curved  horizontal  flame 
front. 


ASYMPTOTIC  ANALYSIS 

While  the  numerical  results  show  us  that  even  a  simple  model  of  a  jet  diffusion  flame  shows 
flame  lift,  our  numerical  results  do  not  show  how  flames  are  stabilized.  In  trying  to  understand 
this,  we  first  briefly  analyze  our  results  in  terms  of  the  simplest  premixed  stabilization  theory, 
which  suggests  that  a  flame  is  stabilized  where  the  cold  flow  velocity  equals  a  premixed  burning 
velocity.  Chung  and  Lee  [20],  for  example,  consider  similarly  simplified  governing  equations  but  for 
an  axisymmetric  jet  geometry  without  a  coflow.  They  obtain  liftoff  heights  by  determining  where 
the  cold  flow  velocity  contour  representing  an  appropriate  premixed  burning  velocity  intersects  the 
stoichiometric  concentration  contour.  We  may  attempt  to  apply  this  methodology  to  the  present 
problem  in  order  to  determine  liftoff  heights.  Since  we  have  assumed  a  one-dimensional  velocity 
field,  the  cold  flow  velocity  in  our  model  is  constant  everywhere  in  the  flow  field.  Then  the  fully 
premixed  stabilization  theory  suggests  that  we  should  find  a  stabilized  flame  for  only  a  single  jet 
velocity  (the  appropriate  premixed  burning  velocity),  with  the  liftoff  height  for  this  velocity  being 
indeterminate  since  the  cold  flow  velocity  is  the  same  everywhere  in  the  flow  field.  Contrary  to 
this,  however,  we  find  stabilized  lifted  flames  for  a  continuum  of  jet  velocities  with  an  apparently 
unique  lift  off  height  for  each  jet  velocity.  This  suggests  that  the  mechanism  is  more  complicated. 
In  an  attempt  to  gain  further  insight  into  the  stabilization  mechanism,  we  therefore  turn  to  theory. 

Specifically,  we  return  to  Eqs.  2-4,  6-8  and  seek  an  asymptotic  solution  in  a  distinguished  limit 
of  large  activation  energy  9  and  large  Damkohler  number  D.  Variables  are  expanded  as 


T  ~ 

To  +  cTi  +  *  *  * 

(19) 

Yf  ~ 

Ypo  +  +  • '  • 

(20) 

Yo  ~ 

Yoo  +  ^Yoi  H - 

(21) 

where  e  =  T‘jl9  is  an  inverse  activation  energy  parameter,  with  Tj  the  flame  temperature.  Solutions 
corresponding  to  a  near  extinction  regime  in  which  fuel  and  oxidant  are  nearly  completely  consumed 
at  the  flame  will  be  sought.  In  such  a  limit  chemical  reaction  is  confined  to  a  thin  zone,  to  be 
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determined  from  the  analysis.  This  “surface”  of  chemical  reaction  will  be  denoted  by  F(x,  y)  =  0. 
To  leading  order,  we  assume  that  the  flame  separates  a  region  of  no  oxidant  (F  >  0)  from  a  region 
where  there  is  no  fuel  (F  <  0).  The  nonlinear  chemical  reaction  terms  in  Eqs.  2-4  will  thus  be 

replaced  by  conditions  at  the  flame  front  determined  by  an  appropriate  analysis  of  the  reaction 
zone  structure. 

The  governing  equations  to  be  solved  at  each  order  are  then 


dT 

d^T 

d^T 

=  0 

^dx 

dx^  ■ 

dy'^ 

(22) 

d'^Yp 

d^Yp 

=  0 

dx 

dx'^ 

dy^ 

(23) 

dYo 

d'^Yo 

d^Yo 

=  0 

dx 

dx^ 

dy^ 

(24) 

which,  again,  are  to  be  solved  separately  for  F  >  0  and  F  <  0,  subject  to  the  boundary  conditions 
Eqs.  6-8,  and  conditions  at  the  flame  front 


Vfo  =  0  for  F  <  0,  Yoo  =  0  for  F  >  0. 


[To]  =  [y>o]  =  [Yoo]  =  0 


(25) 

(26) 
(27) 


[Ti] 

'dTi  dYpi- 

.dn  dn  _ 


-[Xfi]  =  -[Yoi] 
im^dYox] 

’  [dn^  dn 


=  0 


(28) 

(29) 


- ‘5'f  Yoi\f=o-=Sx.  (30) 

Here  jump  conditions  were  obtained  by  integrating  across  the  reaction  zone  [24];  brackets  indicate 

the  jump  m  the  enclosed  quantity  from  F  =  O'  to  F  =  0+,  where  the  normal  derivative  and  unit 
normal  vector  are  given  by 

d  VF 

dn  |VF| 

In  Eq.  30,  Sp  and  Sx  quantify  the  small  amount  of  fuel  and  oxidant  that  leak  through  the  flame 
as  a  result  of  departure  from  complete  combustion  at  0{e).  These  quantities  are  determined  from 
solutions  of  the  flame  front  inner  structure  equation,  to  be  discussed  below.  The  above  conditions 
are  sufficient  to  determine  T,  Yp,  Yq  to  C?(e),  as  well  as  the  flame  position  F(x,  y)  =  0. 

In  general  the  above  system  is  a  free  boundary  problem  which  must  be  solved  numerically 
However,  the  existence  of  the  coupling  functions  Eqs.  11-13  allows  analytical  solutions  to  be  written 
for  the  leading  order  system.  Regardless  of  the  magnitude  of  the  chemical  reaction,  if  at  the  flame 
Jeet  ^el  and  oxidant  to  leading  order  are  completely  consumed,  the  flame  must  be  located  where 
^  ~  0.  Then  from  Eq.  13,  the  flame  position  is  given  by 


.  l  +  Xooil  —  L) 

F{x,  y)  = - ^ - L  +  (1  +  Xoo)S{x,  y)  =  0  (32) 

Resulting  flame  shapes  for  a  similar  flow  configuration  have  been  discussed  previously  in  [14]. 
The  determination  of  the  flame  sheet  now  allows  To,  Ffo,  Too  to  be  calculated  everywhere  in  the 
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combustion  field  from  the  coupling  relations  Eqs.  11.13,  by  substituting  in  the  expansions  Eqs.  19- 
21  and  setting  Ffo  =  0  in  F  <  0  and  n,„  =  0  in  F  >  0  appropriately.  The  adiabatic  flame 
temperature  may  also  be  calculated  from  Eq.  11  or  12, 


ji  _  Too  +  Xoo(l  -f-  Tj) 

f  1  1  ar*  • 


1  +  Xoo 


(33) 


thefl«rp^l^  exploiting  the  use  of  coupling  functions  and  assuming  near  complete  combustion  at 

of  the  frep  be  determined  without  recourse  to  a  fully  numerical  solution 

f  the  free  boundary  problem.  Leading  order  conditions  at  the  flame  and  the  flame  boundary  itself 
are  also  now  known  and  can  be  summarized  as 

atF  =  0-andF  =  0+,  To  =  Tf,  Yfo=^0,  Yqo  =  0  (34) 

where  F{x,y)  and  T/  are  given  by  Eqs.  32  and  33. 

Proceeding  to  0(e),  corrections  to  the  leading  order  solution  are  sought  in  which  departure  from 

reTw  S"llTg”ivT' 


Ti  +  Ypi  —  0,  Ti  +  Yoi  =  0,  Yf\  -  Yoi  =  0 
so  that  jump  conditions  (28)  are  satisfied  automatically,  and 

at  F  =  0-  ,  -Ti  =  Yfi  =  yoi  =  Sf 
at  F  =  0+  ,  -Ti  =  Yfi  =  yoi  =  5;^ . 


(35) 


(36) 

(37) 

rr*”"  ““am  leakage 

through  the  flame.  More  sigiuflcantly  in  the  present  context,  such  an  analysis  will  provide  impor- 
tant  information  regarding  flame  extinction. 

Flame  Front  Inner  Structure 

To  examine  the  reaction  zone  of  the  flame,  coordinates  are  stretched  at  F  =  0  and  inner 
expansions  are  introduced  of  the  form 


F  ~  T/  +  eri  +  ■ 
Yf  ~  €J//  +  • 
To  ~  eyo  + 


(38) 

(39) 

(40) 


^  expansions  into  the  governing  equations  and  performing  an  appropriate  coor- 

in  a  form  firlr™  ""h  T-  h  Ji  temperature  perturbation  can  be  written 

in  a  rorm  first  given  by  Linan  [15]  as 


—  =  (/  _  ^2)  +  tO] 

where  C  is  a  stretched  coordinate  normal  to  the  flame  surface.  Boundary  conditions  are 
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The  conditions  at  the  flame  front  Eq.  30  then  follow  from  the  matching  relations 
=  lim  {<I>  +  Q  =  lim  i4>-Q- 

C->-oo  C-^+oo 

The  reactant  leakage  at  the  flame  sheet  thus  depends  on  two  parameters,  7  and  S,  which  follow 
from  the  coordinate  transformation  and  are  given  by 
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The  parameter  7  is  a  ratio  of  the  excess  heat  transported  to  one  side  of  the  flame  to  the  total 
heat  generated  at  the  flame.  The  limit  of  nearly  complete  combustion  considered  here  implies  that 
-1  <  7  <  1.  When  I7I  >  1  there  is  a  heat  flux  from  the  ambient  to  the  flame,  and,  correspondingly, 
0{1)  rather  than  0{e)  leakage  of  fuel  or  oxidant  occurs  through  the  reaction  zone;  such  a  limit  will 
not  be  discussed  in  this  work.  For  the  particular  configuration  considered  here,  we  find 


_  2(Tj  -  Too)  +  (1  -  Xqo) 

(l  +  Xoo)  ■ 

The  parameter  5  appearing  in  Eq.  41  is  a  reduced  Damkohler  number,  and  thus  a  measure  of  the 
corripleteness  of  the  reaction.  The  limit  <5  — >•  oo  corresponds  to  the  Burke-Schumann  limit  {D  oo) 
of  complete  combustion,  in  which  Sp  =  Sx  =  0.  Numerical  solution  of  Eqs.  41-42  indicates  that, 
for  moderate  values  of  6,  the  system  Eqs.  41-42  has  two  physically  relevant  solutions;  two  different 
quantities  of  the  reactant  leakage  are  possible.  However,  as  5  is  decreased  further,  there  exists  a 
critical  value  6^,  below  which  no  solutions  to  Eqs.  41-42  exist.  This  value  of  5c  depends  solely  on 
7.  Linan  [15]  provides  an  approximation  for  Sc  based  on  numerical  computations  as 


Sc  =  exp(1.0)[(l  -  I7I)  -  (1  -  |7|)2  -f  0.26(1  - \j\)^  +  0.055(1  -  |7|)^] 
This  implies  that  portions  of  the  flame  sheet  where  S  <  Sc,  or 


(47) 


TO]  2  4e^De-^/Tf 

.dn ,  ^  TjSc 


(48) 


are  extinguished. 

From  the  leading  order  solutions  derived  from  Eqs. 11-13,  the  jump  of  the  normal  derivative  of 
the  temperature  across  the  reaction  zone  is  given  by 


[:^J  ((fe) +(^)  («) 

Parenthetically,  we  note  that  this  quantity  is  essentially  the  scalar  dissipation  rate  found  when  a 
mixture  fraction  formulation  is  used  for  the  analysis.  Evaluating  this  quantity  we  And  that  S  varies 
continuously  along  the  flame  sheet,  attaining  its  minimum  value  at  the  flame  base.  Therefore,  for 


7 


moderate  values  of  A  =  Ae^De-^/^f  jT],  extinction  of  the  flame  near  the  flame  base  occurs  below 
X  =  Xc  where  <5(xc)  =  6c. 

When  a  portion  of  the  flame  sheet  is  thus  extinguished,  modification  of  the  leading  order 
solution  is  required.  No  longer  are  there  two  distinct  regions  of  the  flow  field;  in  some  regions  fuel 
and  oxidant  coexist  even  at  leading  order.  Nonetheless,  chemical  reaction  in  these  regions  will  be 
negligible  if  the  temperature  is  sufficiently  low.  To  modify  the  preceding  formulation  to  account 
for  such  a  region,  the  condition  Eq.  25  may  be  relaxed  to 

atE  =  0-,  yFo  =  0,  atF  =  0+,  Too  =  0.  (50) 

where  the  flame  sheet  location  F{x,  y)  =  0  is  now  implied  to  exist  only  where  S  >  6c.  This  relaxation 
of  the  constraint  Eq.  25  does  not  affect  the  validity  of  the  preceding  formulation  if 

^  ^  f  =  0+,  ^  =  0(<).  (51) 

For  D  sufficiently  large  such  that  no  portion  of  the  flame  sheet  is  extinguished,  applying  the 
constraint  Eq.  50  results  in  the  solution  obtained  by  applying  Eq.  25,  with  Eq.  51  being  trivially 
satisfied.  If  the  constraint  51  is  not  met,  the  formulation  must  be  adjusted  and  one  expects  a 
modification  of  the  extinction  criteria  Eq.  48. 

When  solving  the  asymptotic  system  under  the  assumption  51,  for  moderate  values  of  D  one 
must  in  general  still  resort  to  a  fully  numerical  solution  of  the  problem.  In  this  case  the  region  of 
pure  mixing  at  the  base  of  the  flame  will  likely  modify  the  temperature  gradients  at  the  edge  of 
the  flame,  and  thus  the  reduced  Damkohler  number  6,  and  flame  position,  now  defined  as  where 
Yfo  =  Yoo  =  0  and  6  >  6c.  However,  in  a  first  attempt  to  construct  such  solutions,  we  ignore  the 
effect  the  region  of  pure  mixing  has  on  the  flame  downstream,  and  thus  assume  the  constraint  Eq.  51 
remains  satisfied,  and  that  the  extinction  criteria  does  not  need  to  be  modified.  That  is,  conditions 
at  the  flame  are  assumed  to  be  unaffected  by  the  frozen  flow  region;  the  effect  of  the  extinguished 
flame  base  on  the  flow  field  surrounding  the  flame  will  be  taken  into  account.  In  making  this 
assumption,  the  parameters  6  and  7,  as  well  as  the  position  of  the  flame  sheet,  restricted  to  where 
6  >  6c,  remain  as  given  in  Eqs.  45,  46,  32. 

Since  the  solutions  of  the  coupling  functions  Eqs.  11-13  remain  valid  independent  of  the  existence 
or  extinction  of  the  flame  sheet,  any  two  of  the  variables  T,  Vp,  and  Vo  may  be  expressed  as  a 
function  of  the  remaining  one.  Therefore,  again  it  suffices  to  numerically  solve  for  the  temperature. 
The  governing  equation  is  then 

dT  d^T  d^T 

^  dx  dx"^  ~  (52) 

Boundary  conditions  at  y  =  0  and  y  —  L  are 


while  at  a:  =  0, 


(53) 


for  0  <  y  <  1,  mT  - 

The  condition  of  boundedness  as  x 
that  at  X  =  ^maxi 


j^  =  mTr,  forl<y<T,  mT  -  —  =  mToo.  (54) 

-)■  00  will  be  replaced  by  a  condition  on  the  truncated  domain 


(55) 
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(56) 

(57) 


Conditions  at  the  flame  front,  combining  Eqs.  34,  36,  and  37,  are  finally  given  by 

for  (^  >  4  and  F  ^  0^,  T^Tj-  eSp 
for(5><JeandF  =  0+,  T  =  Tj-eSx 

where  F(ar,y),  S  and  7  remain  as  given  in  Eqs.  32,  45,  and  46. 

Solution  of  the  Asymptotic  Model 

The  system  Eqs.  52-57  requires  solution  of  a  steady-state  elliptic  equation.  We  have  solved  the 
system  by  a  simple  five-point  finite-difference  scheme  using  Gauss-Seidel  iteration.  Discretization 
errors  are  thus  of  0{dx^),  where  we  have  taken  dy  =  dx.  Near  the  flame  boundary  F{x^y)  =  0, 
the  discretization  scheme  has  been  modified  to  account  for  situations  in  which  the  flame  does  not 
exactly  fall  on  grid  points;  errors  near  this  boundary  therefore  may  be  slightly  larger. 

Typical  solution  contours  are  shown  in  Figure  7.  Here  Tj  =  Too  =  0.1,  Xqo  =  1.0,  m  =  9.0, 

€  =  0.072  {9  =  5.0),  A  =  1.0,  L  =  10.0,  and  dx  =  dy  =  0.02.  The  initial  condition  for  the  iteration 
was  taken  from  the  Burke-Schumann  solution  for  x  >  Xc  and  the  initial  {x  =  0)  condition  for  x  <  Xc- 
Profiles  show  a  liftoff  height  of  between  one  and  two  jet  diameters,  and  are  in  fairly  good  qualitative 
agreement  with  previous  experimental  and  numerical  studies  [7].  Figure  7  is  qualitatively  similar  to 
Figure  2  obtained  from  the  numerical  solution.  In  the  fully  numerical  solution,  however,  the  flame 
width,  particularly  at  the  base,  appears  broader  than  in  Figure  7,  in  which  the  flame  is  assumed  to 
be  confined  to  a  surface.  However,  this  broader  flame  width  is  consistent  with  a  smaller  effective 
Damkohler  number  at  the  flame  base.  Note  also  that  while  liftoff  heights  are  comparable  in  the 
two  figures,  they  are  obtained  for  largely  disparate  values  of  A  (or  D.) 

Figure  8  shows  the  effect  of  the  mass  flow  rate  on  the  flame  shape.  Here  solutions  of  F(ar,  y)  =  0, 
5>5c  are  shown  for  several  values  of  m,  with  Tj  =  Too  =  0.1,  Xoo  =  2.0  T  =  10.0,  A  =  2.0, 
N  =  200.  As  m  is  increased,  the  liftoff  height  of  the  flame  and  the  overall  flame  height  increase. 
The  length  of  the  flame  (the  distance  from  the  lifted  base  to  the  flame  tip)  increases  as  well. 
In  Figure  9,  flame  shapes  for  the  same  parameters,  but  a  smaller  value  of  the  ambient  oxidant 
concentration  Xpo  =  0.5,  are  shown.  The  liftoff  and  flame  heights  exhibit  the  same  general  trends 
as  in  Figure  8.  However,  for  the  smaller  values  of  Xoo-,  liftoff  heights  are  shorter,  flame  lengths  are 
longer,  and  flame  shapes  are  wider  and  more  rounded. 

In  Figure  10,  liftoff  heights  are  shown  versus  m  for  several  values  of  the  parameter  A.  Here 
Xoo  =  1-0,  with  Tj  —  Too  =  0.1,  L  =  10.0,  and  N  =  200.  As  expected,  liftoff  heights  decrease 
as  A  increases,  consistent  with  the  Burke-Schumann  limit  A  — )■  oc  in  which  flames  are  attached 
for  all  values  of  m.  In  Figure  10  the  flame  height,  which  does  not  vary  with  A  in  the  current 
approximation,  is  also  shown,  and  increases  linearly  with  m.  Flame  height  predictions  from  the 
asymptotic  model  are  compared  with  predictions  from  the  numerical  model  in  Figures  3  and  4. 

The  results  from  the  asymptotic  model  show  some  qualitative  agreement  with  the  numerical 
formulation  results.  Liftoff  heights  predicted  using  the  asymptotically-derived  extinction  criteria 
increase  with  jet  exit  velocity  and  decrease  as  the  Damkohler  number  is  increased.  However,  we 
see  a  linear,  rather  than  faster-than-linear,  increase  in  liftoff  height  with  flow  rate  m.  Also,  for  a 
given  value  of  A,  liftoff  height  predictions  from  the  numerical  formulation  are  much  higher  than 
those  from  the  asymptotic  model.  This  suggests  that  for  a  lifted  flame,  coupling  between  the  frozen 
flow  region  adjacent  to  the  jet  exit  and  the  lifted  flame  base  is  significant.  This  is  recisonable  since 
heat  loss  at  the  flame  edge  due  to  the  adjacent  cold  region  would  reduce  the  effective  Damkohler 
number  there,  thus  resulting  in  an  underprediction  of  lift-off  height.  It  is  reasonable  to  think  that 
if  coupling  were  included  in  the  asymptotic  model,  the  cold  flow  region  near  the  jet  would  tend  to 
expedite  extinction  of  the  adjacent  flame  base,  and  thus  a  faster-than-linear  dependence  of  liftoff 
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height  on  flow  rate  would  result,  along  with  better  quantitative  agreement.  However,  that  the 
asymptotic  model  cannot  account  for  partial  premixing  effects  also  might  be  contributing  to  the 
discrepancy.  A  modification  of  the  analysis  to  account  for  some  degree  of  premixing  would  likely 
modify  the  extinction  criteria  Eq.  48  and  also  might  lead  to  better  agreement  between  the  two 
models.  While  lift-off  height  predictions  of  the  two  models  are  different,  we  note  that  the  flame 
heights  predicted  by  both  models  show  excellent  agreement. 

Because  the  flame  heights  predicted  by  the  two  models  agree,  and  the  general  liftoff  trend  of 
increasing  height  with  increasing  flow  rate  and  decreasing  Damkohler  number  exists  in  both  models, 
we  investigate  further  whether  an  extinction  mechanism  plays  a  role  in  laminar  diffusion  flame  lift 
and  stabilization.  In  the  asymptotic  model,  we  obtain  an  extinction  criterion  that  says  aflame  may 
exist  only  for 

A>i.(l  +  W((|f)  (58) 

Now  the  extinction  parameter  A  involves  primarily  the  Damkohler  number  D  and  the  activation 
energy  6  (recall  A  =  We  therefore  investigate  whether  results  from  the  numerical 

formulation  correlate  with  this  asymptotically-derived  extinction  parameter.  In  Figure  11,  Eqs.  14- 
17  were  solved  to  obtain  curves  of  liftoff  heights  versus  m  for  identical  values  of  A.  Here  the 
values  of  D  and  0  vary  greatly  between  curves,  but  have  been  chosen  such  that  in  all  curves  they 
combine  so  that  A  is  fixed  at  the  value  A  =  40.0.  The  results  show  that  as  the  scaled  activation 
energy  0  is  increased,  the  curves  coalesce.  Since  the  asymptotic  model  was  derived  in  the  limit 
of  large  activation  energy,  one  expects  that  as  this  parameter  is  increased,  solutions  should  begin 
to  converge  to  a  single  solution.  In  the  present  results,  this  means  that  if  results  correlate  with 
respect  to  the  parameter  A,  as  it  is  held  fixed  and  0  is  increased,  the  curves  of  lift-off  height  versus 
m  should  begin  to  coalesce.  Since  Figure  11  shows  that  this  coalescing  appears  to  take  place,  our 
results  suggest  that  the  parameter  A  derived  from  the  asymptotic  analysis  is  indeed  a  parameter 
which  characterizes  flame  lift. 


DISCUSSION 

Experiments  indicate  that  for  low  jet  velocities,  diffusion  flames  are  attached  to  the  jet  nozzle. 
As  the  flow  velocity  is  increased,  such  flames  may  lift  off  of  the  nozzle,  so  that  the  flame  base  is 
located  several  jet  diameters  downstream.  As  the  flow  velocity  is  increased  further,  blow  out  occurs 
and  the  flame  is  extinguished.  Lee  et  al.  [9]  suggest  that,  depending  on  the  fuel  type,  degree  of 
dilution,  and  the  nozzle  diameter,  liftoff  and  blowout  may  both  occur  in  the  laminar  regime,  or 
liftoff  may  occur  in  the  laminar  regime  with  blowout  occurring  after  the  transition  to  turbulence, 
or  liftoff  and  blowout  may  both  occur  only  after  a  transition  to  turbulence.  While  there  have  been 
many  experimental  and  numerical  studies  of  diffusion  flame  liftoff  and  blowout  in  turbulent  regimes, 
fewer  studies  exist  which  investigate  laminar  flame  lift. 

Experiments  involving  laminar  propane  jet  diffusion  flames  by  Lee  et  al.  [9]  suggest  that  liftoff 
heights  increase  as  the  jet  exit  velocity  of  the  propane  is  increased,  at  a  faster-than-linear  rate. 
Flame  heights  increase  linearly  with  an  increase  in  the  nozzle  exit  velocity,  with  the  flame  length 
(distance  from  the  lifted  base  to  flame  tip)  approaching  zero  at  blowout.  When  the  fuel  jet  was 
diluted,  the  liftoff  height  increased  for  a  fixed  value  of  the  jet  velocity.  Lee  and  Chung  [6]  have 
also  investigated  the  effect  of  a  low  velocity  oxidant  coflow  on  a  propane  jet  diffusion  flame.  Such 
a  coflow  has  a  minimal  effect  on  attached  flame  length,  but  results  in  a  larger  liftoff  height.  Liftoff 
and  blowout  of  the  jet  flame  occur  at  lower  fuel  jet  velocities  when  there  is  a  coflow.  Chung  and 
Lee  [20],  [6],  in  a  cold  jet  analysis  similar  to  that  of  Savas  and  Gollahalli  [8],  also  have  derived 
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analytical  expressions  for  liftoff  heights,  based  on  the  assumption  that  a  lifted  flame  at  its  base  is 
fully  premixed.  They  postulate  that  a  lifted  flame  is  stabilized  along  the  stoichiometric  contour 
where  the  (cold)  flow  velocity  equals  the  appropriate  burning  velocity  of  the  premixture.  Lifted 
laminar  diffusion  flames  were  also  observed  in  experiments  and  numerical  simulations  by  Xu  et  al. 
[7].  They  considered  a  circular  jet  geometry  with  an  equal  velocity  oxidant  coflow,  using  detailed 
chemistry  and  complex  transport  models  in  their  numerical  simulation.  Their  results  show  a  lifted 
flame  in  which  profiles  of  the  temperature,  fuel,  and  oxidant  are  qualitatively  similar  to  the  profiles 
shown  in  Figures  2  and  6;  liftoff  height  trends  with  respect  to  system  parameters  were  not  discussed. 

The  results  from  the  numerical  model  show  good  qualitative  agreement  with  such  experimental 
results.  Liftoff  height  trends  shown  in  Figure  3  are  qualitatively  the  same  as  those  described  by 
Lee  et  al.  [9].  Profiles  of  temperature  and  fuel  and  oxidant  mass  fraction  from  the  numerical 
model  are  very  comparable  to  those  obtained  by  Xu  et  al.  [7],  and  profiles  from  the  asymptotic 
model  are  comparable  to  such  results  as  well.  Both  the  current  calculations  and  the  calculations 
of  Xu  et  al.  show  the  existence  of  “wings”  of  fuel  outside  of  the  core  of  the  fuel  jet  at  the  flame 
base,  for  example.  In  the  current  work  these  wings  are  more  pronounced  than  in  [7];  however,  this 
is  likely  attributable  to  the  use  of  a  planar  jet  as  compared  to  the  circular  jet  of  Xu  et  al.,  and 
to  the  exclusion  of  a  radial  velocity  component.  Flame  heights  predicted  by  the  numerical  and 
asymptotic  models  show  excellent  agreement.  That  the  numerical  results  appear  to  correlate  with 
the  extinction  parameter  A  from  the  asymptotic  analysis  suggests  that  diffusion  flame  extinction 
plays  a  role  in  lifted  flame  stabilization. 

It  is  useful  to  again  analyze  the  results  in  light  of  the  flame  stabilization  theory  proposed  by 
Chung  and  Lee  [6].  They  hypothesize  that  flame  stabilization  occurs  along  the  stoichiometric 
contour  where  a  cold  flow  velocity  balances  the  stoichiometric  burning  velocity.  In  the  present 
model,  we  consider  only  a  one-dimensional,  axial  flow  velocity.  As  a  result,  the  velocity  of  a  cold 
flow  is  constant  everywhere  in  the  flow  field.  A  premixed  flame  base  theory  then  suggests  that 
the  model  should  result  in  a  stabilized  lifted  flame  only  if  the  fuel  and  oxidant  source  velocities 
are  equal  to  an  appropriate  premixed  burning  velocity;  in  this  case,  the  flame  stabilization  point 
would  be  indeterminate,  since  at  every  point  along  the  stoichiometric  contour  the  velocity  is  the 
same.  However,  in  the  calculations  we  obtain  unique,  stabilized  lifted  flames  for  multiple  flow  inlet 
velocities.  This  suggests  that  a  propagating,  fully  premixed  flame  base  theory  cannot  fully  explain 
flame  stabilization. 

Our  results  also  should  be  discussed  with  respect  to  partial  premixing  analyses,  [16],  [17],  [18], 
such  as  those  of  edge-flames  and  triple  flames.  In  many  analyses  of  this  type  of  phenomena,  the 
domain  is  considered  to  be  unbounded  and  an  infinite  region  of  partial  premixing  exists.  In  such 
scenarios  it  seems  reasonable  to  equate  a  flame  stabilization  point  with  the  location  at  which  a  flow 
velocity  equals  an  appropriately  calculated  edge-flame,  or  partially  premixed  flame,  propagation 
velocity.  In  the  work  of  Buckmaster  [16]  and  in  that  of  Muller  et  al.  [17],  for  example,  expressions  for 
such  partially  premixed  flame  velocities  are  described  in  which  the  propagation  velocity  is  dependent 
on  the  scalar  dissipation  rate  (which  is  inversely  proportional  to  our  effective  Damkohler  number 
5.)  While  our  results  affirm  the  importance  of  this  parameter,  they  cannot  be  fully  explained 
by  such  a  stabilization  theory  even  when  dependence  of  the  flame  speed  on  scalar  dissipation 
is  taken  into  account.  This  is  in  part  because,  as  the  scalar  dissipation  rate  decreases  and  the 
effective  Damkohler  number  increases  towards  infinity,  such  analyses  reduce  to  a  fully  premixed 
flame  edge  description.  While  such  a  result  may  be  reasonable  and  appropriate  in  the  respective 
models  considered,  clearly  in  the  jet  flame  model  of  the  present  work,  the  infinite  Damkohler  limit 
corresponds  to  a  non-propagating,  attached  Burke-Schumann  flame.  Thus  while  we  cannot  rule  out 
the  possibility  of  partially  premixed  flame  propagation  contributing  to  lifted  flame  stabilization. 
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such  theories  currently  cannot  provide  a  full  description  of  the  phenomena. 

CONCLUSIONS 

A  numerical  solution  of  a  planar  jet  flame  and  a  flame  sheet  model,  in  the  limit  of  large  activation 
energy  and  large  Damkohler  number,  were  used  to  investigate  liftoff  of  laminar  diifusion  flames. 
The  numerical  results  are  qualitatively  consistent  with  experiments  [9]  and  more  complex  numerical 
simulations  [7].  Some  agreement  also  exists  between  experiments  and  the  asymptotic  formulation 
results;  inclusion  of  coupling  and  heat  loss  effects  between  the  frozen  region  below  the  flame  base 
and  the  flame  itself  in  this  model  should  produce  improved  agreement.  Discrepancies  between  the 
two  models  may  be  attributable  to  the  need  to  modify  the  extinction  condition  to  account  for 
flame  edge  effects.  That  the  numerical  model  investigated  here  is  very  simple  and  yet  still  shows 
good  qualitative  agreement  with  experiments  suggests  that  at  least  some  of  the  physical  processes 
responsible  for  flame  lift  are  correctly  accounted  for  in  the  model.  In  particular,  we  note  that  while 
it  is  true  that  ‘^real”  jet  flames  have  two  (or  three)  dimensional  velocity  fields,  it  is  interesting  that 
in  our  numerical  work  we  are  able  to  obtain  lifted  flames  and  triple  flame  structures  with  only  a 
one-dimensional  velocity  field.  Thus  while  the  inclusion  of  multidimensional  eflFects  would  modify 
our  liftoff  results,  it  is  more  significant  that  our  work  clearly  shows  that  flow  redirection  at  the 
flame  base  is  not  necessary  to  obtain  triple  flame  structures  or  flame  stabilization:  The  simplicity 
of  our  numerical  model  in  this  context  produces  insight  which  a  more  complicated  numerical  scheme 
might  be  incapable  of  providing. 

Local  diffusion  flame  extinction  is  shown  to  be  a  plausible  contributing  mechanism  of  flame  lift. 
A  simple  description  of  flame  stabilization  occurring  where  the  flow  velocity  balances  a  premixed 
burning  velocity  is  found  to  be  inadequate.  However,  the  possibility  of  some  degree  of  premixed 
burning  aflFecting  the  flame  structure  and  stabilization  is  not  excluded;  rather,  it  is  shown  that  such 
burning  does  not  appear  to  be  the  sole  controlling  mechanism  of  flame  lift. 

In  turbulent  flames,  Kaplan  et  al.  [19]  have  used  detailed  numerical  simulations  of  the  full 
Navier-Stokes  equations  to  investigate  flame  lift.  Their  results  corroborate  parts  of  both  the  flame 
extinction  and  premixedness  theories,  in  a  turbulent  flame  context.  While  a  turbulent  flow  field  may 
significantly  alter  the  dynamics  of  flame  stabilization,  a  full  solution  of  the  Navier-Stokes  equations 
in  a  laminar  regime  in  which  buoyancy  and  a  multidimensional  velocity  field  are  accounted  for  also 
could  somewhat  modify  the  flame  liftoff  results  obtained  here. 
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Figure  1:  Geometry  of  the  flow  field. 
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Figuie  2:  Solution  contours  of  the  oxidant  and  fuel  mass  fractions  and  the  temperature  from  the 
numerical  model;  m  =  9.0,  Tj  =  Too  =  0.1,  A'oo  =  1.0,  L  =  10,  dx  =  dy  =  0.02,  N  =  200, 
D  =  4  X  10^,  and  0  =  5.0  (A  =  40). 
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flow  rate  m 

Figure  3:  Liftoff  and  flame  heights  versus  the  flow  rate  m  for  several  values  of  D.  Predictions  are 
from  the  numerical  model,  with  Tj  =  Too  =  0.1,  Xoo  =  1-0,  L  =  10.0,  9  =  5.0,  and  N  =  200. 


Damkohler  Number  D 

Figure  4:  Liftoff  and  flame  heights  versus  D  for  several  values  of  the  flow  rate  m.  Predictions  are 
from  the  numerical  model,  with  Tj  =  Too  =  0.1,  Xqo  =  1.0,  L  =  10.0,  6  =  -5.0,  and  N  ~  200. 
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Figure  5:  Temperature  contours  from  the  numerical  model  for  m  =  6,  8,  and  10,  with 
Tj  =  Too  =  0.1,  A'oo  =  1.0,  L  =  10.0,  £>  =  4  x  10^  and  ^  =  5.0  (A  =  40),  and  N  =  200. 
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Figure  6:  Contours  of  the  reaction  rate,  for  m  =  6,  8,  and  10,  with  Tj  =  Too  =  0.1,  Xqo  =  1-0, 
L  =  10.0,  D  =  4  X  10^  and  9  =  5.0  (A  =  40),  and  N  =  200. 
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